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Central Masses and Broad-Line Region Sizes of Active Galactic Nuclei: I. 
Comparing the Photoionization and Reverberation Techniques 

A. Wandel^'2, B.M. Peterson ^, and M.A. Malkan ^ 
ABSTRACT 

The masses and emission-line region sizes of AGNs can be measured by 
"reverberation-mapping" techniques and we use these results to calibrate similar 
determinations made by photoionization models of the AGN line-emitting regions. 
Reverberation mapping uses the light travel-time delayed emission-line response to 
continuum variations to determine the size and kinematics of the emission-line region. 
We compile a sample of 17 Seyfert 1 galaxies and 2 quasars with reliable reverberation 
and spectroscopy data, twice the number available previously. The data provide strong 
evidence that the BLR size (as measured by the lag of the emission-line luminosity 
after changes in the continuum) and the emission-line width measure directly the 
central mass: the virial assumption is tested with long-term UV and optical monitoring 
data on NGC 5548. Two methods are used to estimate the distance of the broad 
emission-line region (BLR) from the ionizing source: the photoionization method 
(which is available for many AGNs but has large intrinsic uncertainties), and the 
reverberation method (which gives very reliable distances, but is available for only a 
few objects). The distance estimate is combined with the velocity dispersion, derived 
from the broad H/3 line width (in the photoionization method) or from the variable 
part (RMS) of the line profile, in the reverberation -RMS method, to estimate the virial 
mass. Comparing the central masses calculated with the reverberation -RMS method 
to those calculated using a photoionization model, we find a highly significant, nearly 
linear correlation. This provides a calibration of the photoionization method on the 
objects with presently available reverberation data, which should enable mass estimates 
for all AGNs with measured H/3 line width. We find that the correlation between the 
masses is significantly better than the correlation between the corresponding BLR 
sizes calculated by the two methods, which further supports the conclusion that both 
methods measure the mass of the central black hole. Comparing the BLR sizes given 
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by the two methods also enables us to estimate the ionizing EUV luminosity Lion 
which is directly unobservable. Typically it is ten times the monochromatic luminosity 
at 5100 A(Lt,). The Eddington ratio for the objects in our sample is in the range 
Lv/LEdd ~ 0.001 - 0.03 and Lion/LEdd « 0.01 - 0.3. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — quasars: 
general — black holes — emission-lines — photoionization 

1. Introduction 

The most fundamental characteristics of the quasar-AGN powerhouse, the central mass 
and structure, are not well understood. The broad emission lines probably provide the best 
probe of these characteristics because the broad-line region (BLR) is the only part of the AGN 

structure where we can measure accurate Doppler motions of gas on appropriate physical scales. 
In particular, assuming the line-emitting matter is gravitationally bound, and hence has a 
near-Keplerian velocity dispersion, it is possible to estimate the virial central mass. 

The virial assumption v oc r~^/^ has been directly tested using preliminary data for only NGC 
5548 (Krolik et al. 1991; Rokaki, Collin-Souffrin, &; Magnan 1993). These works demonstrated 
the case for a Keplerian velocity dispersion in the line-width/time-delay data, for several broad 
lines in two different BLR models. The presently available data support these preliminary results 
and provide strong evidence for Keplerian velocity dispersion in NGC 5548 (Peterson & Wandel 
1999; cf. section 4). Although the virial assumption of Keplerian motion is essential for estimating 
the mass from the BLR observed properties, the estimate may be approximately correct also for 
models in which the line-emitting gas is not bound, such as wind or radiation-pressure induced 
velocities, because the emissivity in a diverging flow decreases rapidly, and most of the emission 
would occure near the base of the flow, when the velocity is still close to the escape velocity, which 
is of the order of the Keplerian velocity (e.g. Murray et al. 1998). 

The main challenges in estimating the virial mass from the emission-line data are to obtain 
a reliable estimate of the size of the BLR, and to relate correctly the line profile to the velocity 
dispersion in that gas. Reliable BLR size measurements are now possible through reverberation 
mapping techniques (Blandford & McKee 1982, recently reviewed by Netzer & Peterson 1997). 
The continuum/cmission-linc cross-correlation function measures the responsivity- weighted radius 
of the BLR (Koratkar &: Gaskell 1991). The spectra can be combined to form the mean spectrum 
and the RMS spectrum, which identifies the variable part of the emission line. This method can 
be used only for a limited number of AGN, because it can be applied only to variable sources and 
requires well-sampled spectra over time scales of months to years. SticIi data arc difficult to obtain 
and exist for fewer than two dozen AGN, mostly low- to moderate-luminosity Seyfert galaxies. 

An alternative method of estimating the virial mass is based on estimating the BLR distance 
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from the central mass using photoionization theory (Netzer 1990), and the emission-hne width 
as an indicator of the velocity dispersion. This method only assumes the line-emitting gas is 
gravitationally bound, but does not depend on the specific BLR geometry or model, so it could 
apply to cloud models as well as disk models (e.g., Rokaki, Boisson, & Collin-Souffrin 1992). 

In its simplest version, the line ratios are used to determine the conditions in the ionized 
line-emitting gas, in particular the density and the ionization parameter U (the number of ionizing 
photons per electron) 

U = ^ = ^ , (1) 

where 

Q = f ^dv 

JlRyd IIV 

is the number of ionizing photons, Ljo„ is the ionizing luminosity and E = Li^jQ is the average 
energy of an ionizing photon. With an estimate of the luminosity of ionizing spectrum, (Mathews 
h Ferland 1987; Bechtold et al. 1987; Zheng et al. 1997; Laor et al. 1998) the eq. |l| can be 
inverted to give the distance of the ionized gas from the continuum source. Wandel (1997) has 
demonstrated that the BLR sizes derived even by the basic photoionization method are in good 
agreement with the reverberation sizes. He also showed that using soft X-ray data to obtain a 
better estimate of the ionizing flux can signiflcantly improve the correlation, quantifying the error 
introduced by the uncertainty in the ionization parameter \J and the density ng. 

A recent compilation of the available objects with reverberation data (Kaspi et al. 1997) 
shows that the BLR size scales roughly as the square root of the luminosity, which agrees with 
the intuitive prediction of the basic photoionization method, assuming that BLR in different AGN 
have similar values (or a narrow, luminosity-independent distribution) of the ionization parameter 
and density. 

In contrast to the reverberation method, the photoionization method is indirect and the BLR 
sizes derived are subject to many intrinsic uncertainties. The line emission is probably extended 
in radius, and includes emission from gas with a range of densities, so that single-zone models 
oversimplify the real situation. While simple photoionization models may not be trusted for 
estimating BLR parameters, it should be possible to calibrate the photoionization method using 
reverberation results (Wandel 1998). 

Since the photoionization method can be used for any AGN, even with low-resolution 
spectroscopy (and even with a single observation), it has enormous potential. If the method can 
be calibrated, as we attempt to in this work, it could give virial masses of a large number of AGN, 
virtually all of those having reasonably well-established H/3 line widths and UV line ratios. A 
growing number of high-redshift quasars are also now being observed in the H/3 line with good 
infrared spectroscopy, which we will also use (Mcintosh et al. 1998). Such a calibration will provide 
us with a powerful tool to determine reliable BLR sizes and virial masses of a large number of 
AGN, and improve our understanding of the MjL relation and the structure of the BLR. 
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In the present work we investigate whether or not such a cahbration is feasible, using ah the 
presently available reverberation and emission-line variability data. 

An important feature of this paper is that we present a sample of AGN that have been 
analyzed in a consistent fashion, providing reliable estimates (to within basically geometrical 
factors of order unity) of AGN black hole masses. The method of estimating virial masses is 
superior to what has been done previously, as we use the scale length from reverberation techniques 
and only the variable part of the emission line to obtain a Doppler width by measuring the line 
width in the rms rather than the mean spectrum. The sample analyzed here is about twice as 
large as any previously discussed sample. 

In the next section we describe the photoionization method. Section 3 describes the 
reverberation - RMS method and data. In section 4, we consider the multi-year observations of 
NGC 5548. Sections 5-6 give the correlation results for the central masses and the BLR sizes, 
section 7 discusses the ionizing luminosity and Eddington ratio estimates, and section 8 discuses our 
result and the uncertainties associated with the sample selection, with the reverberation methods 
and of the photoionization method. 



2. The Basic Photoionization Method 

As a preliminary test we use the calibration of the ionization-parameter method in its most 
basic form — the BLR modeled by a single thin shell, with an average value of UUf, assumed for 
all objects. Combining the definition of the ionization parameter with the assumption of virialized 
velocity dispersion gives the simplest version of the photoionization virial mass estimate: 

M,,.!^ = k{^Ym^ (2) 

where K = fkG~^{A'Kc)^^/'^ and fk is the factor relating the effective velocity dispersion to the 
projected radial velocity deduced from the emission-line Doppler broadening, and Rph is the radius 
derived from eq. |l]. Relating the velocity dispersion to the FWHM we assume < v"^ >= |wfwhm) 
so A = 3/4 (Netzer 1990). 

Since we intend to find the normalization factor with respect to the reverberation-rms direct 
method, we do not have to worry about the precise factors, but keeping the constants can give an 
estimate of the effective values of the physical parameters. The size of the BLR is then derived 
from eq. (^. As we do not know the ionizing luminosity, we need to estimate it. In this work, as 
a zero-order estimate we assume it is proportional to the visible luminosity, and try to determine 
the sample average of the proportionality constant by the calibration. A priori it is reasonable 
that this assumption would contribute to the errors, because individual AGN may have different 
Lion/Ly ratios. However, as we see below, apparently this is not the case, at least for the sample 
at hand. More elaborate estimates of the ionizing continuum may be attempted by using near UV 
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(Zheng et al. 1997) or soft X-ray (Wandel 1997; Laor 1998) data, but it is not clear how much these 
improve the estimate of Lion-, because we do not know the shape of the EUV, and these data are 
also not available for many AGN; our intention is to develop an "affordable" method, applicable to 
as many AGN as possible. Finally, the spectral correction enters as a square root, which reduces 
the eventual errors contributed to the mass estimate. Relating the ionizing luminosity to some 
measured luminosity L eq. (||) gives 

where -L44 = L/ 10^^ erg s~^ and Jl = Lion/L is the factor relating the observed luminosity L to 
the ionizing luminosity, Ei = £'/lRyd which gives 

Mp^, « (2.8 lO^Mo)/ 1^-^^%,^, (4) 

where / = fkf]!^ Ei~^^'^ and V2, = UpwHM/lO^km s"^ 

As a first approximation, we assume the ionizing luminosity is proportional to the visual 
luminosity Ly ~ vLy(b\QQ K). For our purpose this is equivalent to assuming that the ionizing 
luminosity is proportional to uLy. 

As discussed in section 7 below, comparing the reverberation and photoionization BLR sizes 
can be used to estimate fi- 



3. Reverberation-RMS Masses 
3.1. The Method 

The reverberation BLR sizes and the masses are derived by Peterson et al. (1998a)0. The size 
of the BLR is measured by observing the light travel-time delayed response of the emission line to 
continuum variations. The relationship between the light curve L{t) and the emission-lineintensity 
I{t) is assumed to be 

m = J ^{T)L{t - T)dT, 

where ^'(r) is the transfer function (Blandford & McKee 1982). The transfer function depends 
on the BLR geometry, the viewing angle and the line emissivity. The cross-correlation of the 
continuum and emission-line light curves is given by 

CC{t) = f ^{t')AC{t' - T)dT', 



^On account of several typesetting errors in the original paper, the values for masses are superseded by the masses 
given in Table 2 below. 
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where AC is the continuum autocorrelation function (Penston 1991). It can be shown (Koratkar 
&; Gaskell 1991) that the centroid of the cross-correlation function , Tcent, gives the size CTcent 
associated with the emissivity-weighted mean radius for the emission region of a particular 
emission-line. 



The root-mean-square (RMS) spectrum defined as 



ct(A) 



TV 



(Ar_i)-i^(i^(A)-F(A))2 



1/2 



(5) 



where N is the number of spectra and F is the average spectrum. The RMS spectrum measures 
the variations about the mean, automatically excluding constant features such as narrow emission 
lines, Galactic absorption and constant continuum and broad line features. The RMS H/3 
emission-line profile gives the velocity dispersion in the variable part of the gas, the one which is 
used to calculate the BLR size. 

The line width and the BLR size may be combined to yield the virial "reverberation" mass 
estimate (equivalent to eq. ^ 

M.,. . (1.45 X lO^Mo) (i-J^^^) Vrms/, (6) 
where Vrms,3 = Vp^iiMirms) / km s~^. 



3.2. The Data 

Our sample consists of 19 objects, 17 Seyfert 1 galaxies and 2 quasars. This is the most 
complete compilation to date of AGNs for which reverberation data are available. The sample will 
be described in a separate paper, and in this contribution we consider only ground-based optical 
data. These data are drawn from the following sources: 

1. The Ohio State AGN Monitoring Program. A program of approximately weekly 
observations of the H/3 spectral region in bright Seyfert galaxies was carried out with the 
Ohio State CCD spectrograph on the 1.8-m Perkins Telescope at the Lowell Observatory 
from 1988 to 1997. Data on Akn 120, 3C 120, Mrk 79, Mrk 110, Mrk 335, Mrk 590, and 
Mrk 817 used here have been published by Peterson et al. (1998a). The data used here on 
NGC 4051 are preliminary results based on OSU data only, and these are currently being 
combined with data from other sources. 

2. The International AGN Watch. Over the last decade, this consortium has carried out 
multiwavelength monitoring programs on a number of AGNs, and the data are publicly 
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available]^ (Alloin et al. 1994; Peterson 1999). In this study, we employ the AGN Watch 
optical data on 3C 390.3 (Dietrich et al. 1998), Fairall 9 (Santos-Lleo et al. 1997), Mrk SOSQ 
(Carone et al. 1996), NGC 3783 (Stirpe et al. 1994), NGC 5548 (Peterson et al. 1991, 1992, 
1994, 1999; Korista et al. 1995), NGC 4151 (Kaspi et al. 1996a), and NGC 7469 (Collier et 
al. 1998). 

3. The CTIO AGN Monitoring Program. In connection with the AGN Watch program 
on NGC 3783 (Stirpe et al. 1994) in 1992, observations of NGC 3227 (Winge et al. 1995) and 
IC 4329A (Winge et al. 1996) were also made. For the sake of greater internal consistency, 
we reanalyzed these spectra in the same fashion as the Ohio State monitoring data. 

4. The Wise Observatory QSO Monitoring Program. A spectroscopic monitoring 
program on low-redshift quasars is being carried out at the Wise Observatory. Thus far, 
preliminary results on PC 0804+762 and PC 0953+414 have been published (Kaspi et al. 
1996b), and the authors were kind enough to make the original spectra available to us for 
this analysis. 

For each galaxy, we computed the optical continuum~H/3 cross-correlation using the 
interpolation method of Gaskell & Sparke (1986) as implemented by White & Peterson (1994), and 
determined the centroid Tcent of the cross-correlation function using all points within 80% of the 
peak value rmax- Uncertainties in Tcent were evaluated using the model-independent Monte-Carlo 
FR/RSS technique of Peterson et al. (1998b). Spectra were combined to form average and 
root-mean-square (rms) spectra, and the full-width at half maximum of the H/3 feature in the 
average and rms spectra was measured, transformed to the rest frame, and is quoted as a Doppler 
width in kilometers per second. 

A summary of the basic data used in this investigation is given in Table 1. The name of the 
galaxy appears in column (1). Column (2) gives the monochromatic luminosity ALa(5100A) (ergs 
s~^) of the galaxy at a rest wavelength of 5100 A. This is the mean luminosity observed during 
the monitoring campaign from which the data were drawn. In each case, this has been corrected 
for Galactic extinction using the values in the NED database]^. Columns (3) and (4) give the 
full-width at half-maximum of the H/3 emission line in the mean and rms spectra, respectively. 
Column (5) gives the time lag Tcent for the response of the H/3 emission line, and the associated 
uncertainties, which are not symmetric about the measured value. 



^International AGN Watch can be obtained on the World-Wide Web at URL [ittp://www.astronomy.ohio^ 
5tate.edu/~agnwatch / 

^Supplemented by additional data from Peterson et al. (1998a). 

^ The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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We expect the line width in the rms spectra to be larger than the line width in the average 
spectrum because in the former the narrow-line component is eliminated. However there are 
objects for which the FWHM of lines is is narrower in the rms spectrum than in the average 
spectrum. One possible interpretation is that the emission-line gas with the highest velocity 
dispersion and nearest to the central source, which contributes to the far wings does not vary 
much, perhaps because it arises in a medium that is optically thin to ionizing radiation (e.g.. 
Shields, Ferland, & Peterson 1995). 

Note that for NGC 7469 we give the lag with respect to the UV continuum rather than relative 
to the continuum at 5100 A. The reason for this is that in NGC 7469, the optical continuum 
variations follow those in the ultraviolet continuum by about one day. This is the only AGN 
in which this effect has been measured (Wanders et al. 1997; Collier et al. 1998; Peterson et al. 
1998b). 

4. NGC 5548: Evidence for the Keplerian Assumption 

The basic assumption of this model is that the line-emitting gas is gravitationally bound, 
and its velocity dispersion is approximately Keplerian. If this is the case, we would expect that 
during variations of the ionizing luminosity, the "virial mass", given by M ~ G~^v'^r remains 
constant. While the radius given by the luminosity scaling of the photoionization method (eq. ^ 
is uncertain, the reverberation radius and the associated velocity derived from the variable part of 
the line (FWHM-RMS) provide a direct measure of the virial mass. 

The virial assumption can be tested for individual objects. It is known that the BLR is 
stratified, with different emission lines different distances from the central source. Typically 
different lines also have very different line widths. If the motions in the BLR are dominated by the 
gravity of the central mass, the Keplerian relation (P) applied to different lines in the same objects 
must yield the same value for the central mass. Similarly, if the continuum and emission-lines vary 
over long time scales (the short time scale variability is applied for deriving the reverberation BLR 
size), the virial mass should remain constant. In other words, above test states that the lags 
{R) of different lines and at different epochs should be anticorrelated with the rms line width as 
7? oc 

We have tested the virial assumption with the multiyear multiwavelength data on NGC 
5548 (Peterson & Wandel 1999). The emission-line lags r range from ~ 2 to ~ 30 days, and 
are anticorrelated with the line widths FHWM, consistent with a constant value of the product 
r X (FWHM)'^, as expected if the virial assumption is correct. There is significant scatter around 
the best virial fit to the data, however; we suspect that this is due in part to the fact that both 
the lag and the line width are dynamic quantities that are changing in response to changes in the 
continuum flux. Moreover the H/3 results for two years (1993 and 1995) had to be excluded from 
this analysis because the rms profiles were strongly double peaked and therefore did not yield an 
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unambiguous line width. Again this suggests to us that while the virial interpretation is correct to 
first order, it is not the whole story. 

Figure 1 shows the virial mass calculated from the reverberation radius and the line width 

in the rms spectrum as a function of the FWHM-rms for each line in each year. This way of 
presenting the data enables us to demonstrate that the virial mass of all lines is consistent with a 
single value. 

Fitting an M=const. to the points in the M vs. FWHM plane is equivalent to fitting a line 
with a slope of -2 (lagoc(FWHM)^^) in the lag vs. FWHM (log- log) plane (if the two observables 
were uncorrelated, the best fit would be M oc {FW H M)"^) . 

We find an average virial mass (weighted by the uncertainties in the mass) of (6.1 ±2) x lO^M©. 
The weighted mean has been calculated as < x >= S(xi/c7?)/S(l/c7?) where (jj is the uncertainty 
in Xi in the appropriate direction (up or down) with respect to the mean. It is easy to see that 
all the points in the figure are consistent with the average value within ~ Icr. While the annual 
averages of line width and lag in the NGC 5548 data set vary by factors of 2-3 and more than 
an order of magnitude, respectively, they are anticorrelated in a manner close to the Keplerian 
relation above, so that the virial mass remains approximately constant (within the uncertainties). 

We have calculated the goodness of the M=const assumption in three schemes: 

• An M = M fit in the M vs. FWHM plane 

• A y = —2x -I- a fit in the log(lag) vs. log(FWHM) plane (that is, x-y plane) 

• The best fit y = 6a; -|- a in the log(lag) vs. log(FWHM) plane 

where M is the uncertainty- weighted average reverberation mass, x =log(FWHM) and y=log(r). 
We have performed these fits for five data sets: 

• All data points 

• All lines and H/3 represented by the weighted average of the six years (log(FWHM)=3.71, 
log(lag)=1.25) 

• All lines from the years 1989 and 1993 (when available), i.e. all lines and H/3 from 1989 

• Only H/3 for the six years, lags from the centroid of the CC function 

• Only H/3 for the six years, lags from the peak of CC (see below) 

The last two sets demonstrate the constancy of the virial mass for time variability of a single line, 
but evidently the small number of points does not give very significant results. The last set may 
test the sensitivity of the reverberation scheme to the geometry (see below) . 
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For all sets we have calculated the reduced in the three schemes (table ^) The best fit 
and the in the (x,y) plane are calculated with weighting by lag uncertainties (Bevington pp. 
118 and 188). The errors in the x — y plane were estimated from the (logarithmic, directional) 
standard deviation for each point i, 

cTj = log(l + o-+/Tj) if yi<bxi + ai 

ai = - log(l + CF~ /Ti) if yi > hxi + Oj 
where CTi^ and are the + and - uncertainties in the lag r. 

Note that the M=const fits have lower values of than in the FWHM-lag fits (in the x — y 
plane). Although the two are equivalent, the latter uses only the uncertainties in the lag, while 
the former uses the uncertainties in M, which combine the lag and FWHM uncertainties, and is 
therefore more appropriate. 

The best fits for the three sets with all lines are consistent with the virial assumption (slope of 
-2 in the x-y plane). On the other hand, the H/3 fits are not consistent - this is due to their small 
dynamical range, which is even reduced by the weighting, as the point with the largest FWHM 
also has a much larger uncertainty in the lag than the other points. Also the dynamical nature 
of the lag and line width, which depend on the constantly changing continuum flux is probably 
contributing to the scatter. Given these factors and the small number of points, this result is 
probably not very significant. Note however that the tests give a reasonable probability to the 
M=contant assumption even for the H/3 sets, and certainly cannot exclude it. 

For the H/3 annual data we have tried also a different reverberation scheme, using the lags of 
the peak (rather than the centroid ) of the line-continuum CC function. This may give an estimate 
of how sensitive these results are with respect to the uncertainties of the reverberation technique 
and its dependence on the BLR geometry (see section 8.2). 

The (uncertainty-weighted) mean of the annual virial masses derived from the H/3 line annual 
averages is 

< Mrev-c >= (6.6 ± 1.6) X l(fMQ H/? annual, centroid lag 

and 

< Mrev-p >= (6.6 ± 0.5) X W'^Mq H/? annual, peak lag. 

In addition to being mutually consistent, these results are in good agreement with the value 
derived by using all the available emission- lines and years (altogether 14 data points), 

< Mrev >= (6.1 ± 0.4) X 1O^M0 different hues, centroid lag. 
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5. The Mass Correlation 

The basic procedure is quite straight forward: use rms spectra of the objects with available 
reverberation analyses, and compare the virial mass calculated this way to the value calculated 
with the ionization parameter method, using canonical AGN values for Uue and Q/L rather than 
specific values for each object. 

Table 2 compares the BLR sizes and the virial masses derived by the two methods. 

We now try to find the correlation and the functional relation between the masses determined 
by the two methods for the objects with reverberation data. We expect the relation to be linear, 
Mrev oc Mph, but we are also testing for a different power-law dependence. Since the parameters 
that we have approximated appear in the square root, the scatter introduced through these 
approximation may be small. Since we find that the two mass estimates are well correlated, the 
correlation may be used to derive a functional dependence — to a first approximation we assume 
a linear correlation, which gives the normalization coefficient in eq. (^). Figure 2 shows the virial 
masses calculated by the photoionization (eq. |^ and the reverberation-rms methods. 

Considering the unweighted data points gives log(Mj.et,) = 0.82 log(Mp^) + 1.53 with a 
correlation coefficient of 0.83. 

5.1. Weighting by Measurement Uncertainties 

Since the uncertainties in Mrev are very different for different objects, it is important to weigh 
the points according to their uncertainties. Taking into account the directed uncertainties in the 
Mrev the correlation coefficient becomes 0.89, and the best fit 

log{Mrev/MQ) = b log{Mph/MQ) + a, 

with b = 0.93 ± 0.07 and a = 0.70 ± 0.53. 

For NGC 5548 masses we have used the weighted averages of the six years. In the weighted 
fit we have increased the formal uncertainty of the Mrk 590 and NGC 5548 reverberation mass 
estimates by a factor of ~ 3, as they were significantly smaller then the other objects, and gave an 
excessive weight to these points. 



^Directed in the sense of using the appropriate error bar (upper or lower) of each point with respect to the fit, 
and repeating the fitting until it converged. 
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5.2. Weighting by Intrinsic Uncertainties 

We have used the standard procedure of weighting data with uncertainties (Bevington 
1969, p. 118) which takes into account only uncertainties in one (the dependent) variable. This 
weighting method may give too much weight to points with a very small uncertainty in Mrev, 
since the uncertainty in Mp^ is always large, as it has to account for the intrinsic uncertainties in 
the photoionization method, discussed above. One could apply a "mixed weighting covariance" 
method, where each data point is weighted considering its x and y uncertainties (a detailed 
discussion of this subject can be found in Feigelson and Babu, 1992). However, also the 
reverberation method has intrinsic uncertainties, as discussed in section 8 below. The errors 
quoted for Mrev represent only the formal measurement errors, given by 



AM 



1/2 

(7) 



Therefore also M^ev has an additional intrinsic uncertainty (possibly larger than the measuring 
errors). Taking into account significant intrinsic uncertainties in both, M^ev and Mp/^ has the 
effect of giving all points a similar weight, i.e. an unweighted fit. 



5.3. Interpreting the Mrev — Mph Relationship 

Our most important result is the very good correlation, consistent with a linear relation, and 
the surprisingly low scatter. We note that the outlying objects in our sample are less than a factor 
of two away from the diagonal {Mp^ = Mrev) line. 

The slope (6) in the logarithmic fit deviating from unity may account for intrinsic residual 
dependence of the photoionization mass on on parameters such as luminosity or line width, (for 
example, if one of the parameters in eq. (^, such as U or Lion/Ly were dependent on L or on M). 
The fact that we find the slope to be consistent with unity implies that there is no evidence for 
such an effect in our sample. 

The constant (a) in the logarithmic fit gives an estimate of the calibration coefficient, if b is 
close to unity. A more direct determination the calibration coefficient is obtained by imposing 
6 = 1 and solving for the best value of a, which we denote by ai. If 6 is close to unity, this value is 
given by 

Sm =<y> - <x>, (8) 

where < x >= x and < y >= y are the averages (weighted, Bevington p. 88, or unweighted, 
as discussed above). For our sample we find < log{Mrev/MQ) >= 7.49 and 7.64 while 
< log{Mph/MQ) >= 7.29 and 7.48, for the unweighted and weighted averages, respectively, so that 
6m = 0.20 and 0.16, respectively. 
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The difference 6m is actually the empirical factor which best calibrates the 
photoionization method for determining the virial masses, so for the present sample and 
our choice of parameters we find for the actual coefficient in eq. (^: 

/ = 10''" w 1.6 and 1.45, 

respectively. 

6. The Radius Correlation 

The errors in the photoionization virial mass determination (treated in more detail in 
the discussion section) consist of two main ingredients: the uncertainties introduced by the 
photoionization method of estimating the BLR size and the uncertainty in matching the size 
derived to the velocity dispersion implied by the line-width (and which measure of the line profile 
best describes the relevant velocity dispersion). While these two problems are automatically solved 
in the reverberation - rms method, it important to understand their contribution to the errors in 
the photoionization method. 

By comparing the BLR size estimates of the photoionization method with the reverberation 
radii we can estimate which part of the scatter in the Mph — M^ev correlation (Fig. 2) is due to 
the uncertainties in the photoionization model approximation and the assumptions of average 
ionization parameter, density, and Lion ~ Ly- 

The empirical result rrev oc L^/^ (Kaspi et al. 1997) would predict a linear relation, 
\og{rrev) = l-01og(rp/,) + C, as the photoionization method in its simplest form implies oc L^/^. 

The data and their interpretation are shown in figures 3a and 3b, respectively: Figure 3a 

shows the lag vs. the monochromatic luminosity at 5100A. The logarithmic linear (unweighted) fit 

is loglag(days) = (0.46 ± 0.1)logL44 + (1.50 ± 0.11) where L44 = ALA(5100l)/10^^erg s'^. Within 

the uncertainties this is consistent with the result of Kaspi et al. (1997). Forcing an i? oc L^/^ slope 

1/2 

gives an almost identical result, Rhp = 33L44 light days, or, expressed in terms of the luminosity 
integrated over the 0.1-l//m band rather than the monochromatic luminosity (assiming L;^ oc A or 
oc u-^) 

Rhp = 15L4{^(o -^^-^^j^^hght days. 

Figure 3b shows the lag vs. the BLR size estimated by the photoionization method for our 
sample. The logarithmic linear fit is lag(days) = 0.92r(lt-days) + 0.45. The correlation (r = 0.70) 
is significantly less good than the correlation of the masses. This is caused by a few objects 
with large uncertainties being well off the diagonal. When the points are weighted by the lag 
uncertainties, the slope becomes more deviant from unity — b = 1.3 it 0.12 — because several 
off-diagonal objects with small uncertainties dominate the fit. As discussed above, the unweighted 
fit may be more realistic as it approximately gives the result of adding the intrinsic uncertainties. 
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Efectively it softens the influence of the objects with very small lag error bars in the weighting, 
lowering the slope back to be consistent with unity. 

It is interesting to determine the calibration coefficient also for the BLR size, derived from 
the photoionization method. Considering figure 3b, the weighted averages are < log(lag) >= 1.22 
and 1.31 and < log{rph) >= 0.81 and 0.94 for the unweighted and weighted methods, respectively. 
The calibration factor =< log(lag) > — < \og{rph) > is then 0.39 and 0.37, so that the radius 
calibration factor is 

fR ^ lOfj = 2.4 (9) 

for both weighting methods. We see that the correlation between the BLR sizes given by 
the reverberation and photoionization methods is significantly weaker than the corresponding 
correlation between the masses and the scatter is larger. 



7. Estimating Lion and L/M 

7.1. The Ionizing Luminosity 

The coefficient in eq. @ is / = jfj/'^Ei ^^"^ . We may reverse this equation and use the 
value we have found for / in order to calculate the effective ratio /l between the visual luminosity 
ALa(5100A) and the ionizing luminosity: 

However, the empirical factor / actually accounts for two components: the luminosity factor, //^, 
and the kinematic factor, f^, which we have taken to be | in both methods, but it may contain 
also a possible factor contributed by the mapping from the FWHM in the photoionization method 
to the FWHM-RMS in the reverberation method. 



A more direct determination of /l can be obtained from the calibration factor of the BLR 

1/2 1/2 

sizes: since Rph oc L^j,^ while the true emission-line radius measured by the lag is Rrev oc L-^,^ , we 
have Lion/Lobs {Rrev / Rphf ■ In our notation, eq. (^) gives 



which for our choice, U nio = 1 becomes 

Jl = lO^-^^^i = 6^1 

. For reasonable ionizing spectra (Wandel 1997, Bechtold et al. 1987) E ~ 1.2 — 2, so that that for 
our sample the ionizing continuum may be related to the monochromatic luminosity at 5100 A by 

Uon ~ 10?7n7-i2 ALa(5100A). 
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For example, we can compare this empirical estimate of the average photon energy to some 
published model estimates. Integrating the "Medium" AGN spectrum described in Bechtold et 
al. (1987) gives a total bolometric luminosity of 1.4 times the value of i^L^ at 1 Rydberg, which 
equals 7.7 times the value of vL^ at 5100 A. The exact values are relatively insensitive to the 
upper limit to the energy integral (the choice of hard-X-ray cutoff). About half of this total power 
is emitted by the "power-law" component longward of 1500 A. For Bechtold et al.'s "Hard" AGN 
continuum energy distribution, a somewhat smaller fraction of the total bolometric power emerges 
at wavelengths longer than 1 Ryd, and the total bolometric power is 13 times the value of pLy at 
5100 A. Within the uncertainties on average BLR parameters, either of these models is consistent 
with the average AGN continuum spectrum inferred from our comparison of photoionization and 
reverberation estimates of BLR size. 



7.2. The Eddington Ratio 

We may use the high quality virial mass estimates obtained by the reverberation -RMS 
method for this sample, to estimate a mean L/M relation. Combining this with the above estimate 
of the ionizing luminosity (which presumably is comparable to the total bolometric luminosity, 
as a large part of the AGN energy is expected to be radiated in the UV) we may estimate the 
Eddington ratio. 

Figure 4 shows the reverberation - RMS virial mass vs. the 5100 A luminosity for our sample. 

An unweighted linear fit gives log(M/M0) = 0.54 log(L44) -|- 7.94, with a correlation coefficient 
r = 0.68. However, here the fit needs to be weighted only by the mass uncertainties, as those of 
the luminosity are relatively small. A weighted fit yields 

log(M/M0) = (0.77 ± 0.07)log(L44) + (7.92 ± 0.04). (11) 

We can calculate the Eddington ratio from the difference of the weighted averages, 

Sml =< log(M/M0) > - < log(L44) >= 8.0. (12) 

Combining this with the definition of the Eddington luminosity LEdd = 1-3 x 10^^ {M/Mq) erg s~^ 
we have for our sample 

ALa(5100A) «i O.OlLEdd (13) 

and 

Lion/LEdd^O.lUniQ. (14) 
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7.3. Ionizing luminosities of individual objects 

We may apply eq. |l^, (or inverse eq. (3)) to estimate the ionizing luminosity of individual 
objects in our sample. This gives 

Lion-EiUnioXLxi5WOA)(^^^ = 6 x lO^^^iC/nio (^^^^) erg s^\ (15) 

where r is the lag (the centroid of the cross-correlation function) . Since the ionizing luminosity is 
likely to be a large fraction of the bolometric luminosity, and the virial reverberation -rms mass is 
our most reliable estimate of the mass of the central black hole, these values give a direct measure 
of the true Eddington ratio 

L/L^dd ~ Lion/Mj-ev 

The last two columns in table 2 show the ionizing luminosity and total Eddington ratio of the 
objects in our sample. 

Fig. 5 shows the virial (reverberation) mass vs. the ionizing luminosity for the objects in our 
sample. We note that most objects in our sample have L/LEdd ~ 0.01 — 0.3. We also note that 
there is a strong trend of the Eddington ratio to increase with luminosity: an unweighted fit to 
the Eddington ratio vs. luminosity gives 

L/LEdd = O.OSL^f , 

where L44 is the luminosity at 5100A. However, since Ljon oc (eq. p!5| ) and M oc r x [FWHM)"^ , 
this correlation may be a mere reflection of the correlation of the lag with luminosity. 



8. Discussion 
8.1. Results 

From the present sample we find that the virial masses of the 19 AGN in our sample estimated 
by the photoionization method are very well correlated with the more direct estimates from 
the reverberation - RMS method. We find the relation is consistent with a linear one, that is, 
Mrev cx Mph. For an average value of Uue = 10^° and assuming Lion = /lAFx(5100A) we find a 
calibration factor of / ~ 1.5 where Mrev = fMph, which corresponds to Jl ~ 10. 

The fact that the photoionization masses are better correlated with the reverberation - RMS 
masses than the corresponding BLR size with the reverberation distances contradicts the intuitive 
prediction that the scatter in the masses should be larger, due to the additional uncertainty 
in translating line profiles to the appropriate velocity dispersion. This may indicate that there 
is an internal compensation mechanism, which reduces scatter in the photoionization mass 
determination. This may be related to the difference between the reverberation and the 
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photoionization methods: the former measures the size of the variable component of the BLR, 
while the latter measures the emissivity-weighted size, and the two may differ. However, when 
applied to the virial mass, both methods measure the same mass. 

This conclusion is strongly supported by the variability of NGC 5548 (Fig. 1), where the 
variations in the lag of the line width are significantly larger than the variations in the combination 
lagxt'p^jjj^. (The latter is actually consistent with remaining constant, within the uncertainties). 

Also, examining the deviations of individual objects from the linear relation in the mass 
and size plots, we note that several objects which deviate significantly in Fig. 4 (the BLR size 
correlation) have smaller deviations in the mass correlation (Fig. 3). 

We have also estimated the ionizing luminosity, the mass-luminosity relation and the 
Eddington ratio for our sample. Applying the calibrated photoionization method to many AGN 
may reveal a trend for these important parameters in various groups of AGN. 

Obtaining reliable black hole masses for many AGN may also provide important insight for 
the relation between AGNs and their host galaxies, as has already been demonstrated for the 
present sample (Wandel 1999). 

8.2. Comparison with Accretion-Disk Models 

AGN masses may be estimated also by fitting the continuum optical- UV spectrum by a 
thin accretion disk model spectrum (e.g. Malkan 1990; Wandel and Petrosian 1988; Laor 1990). 

We have compared the reverberation and photoionization mass estimates with those obtained 
by standard thin-accretion disk fitting of the continuum energy distribution. Specifically, more 
than half of the calibrating Seyfert I's have been fitted in Sun and Malkan (1989), Malkan 
(1990), and Alloin et al. (1995). The black hole masses estimated from continuum fitting are 
generally correlated with those derived from the emission lines. We find that in general the 
black hole mass obtained from the accretion disk method tends to be larger than the mass 
found from the reverberation data. In order for the normalizations to be consistent with the 
reverberation masses the parameters of the disk models have to be chosen to give the smallest 
possible black hole masses. Specifically, the disk fits must all assume non-rotating (Schwarzschild) 
black holes, and zero inclination (face on) to give a normalization consistent with the emission 
line estimates. For most of the objects in our sample, inclined disks, or rotating (Kerr) black 
holes would yield larger continuum- fitting black hole masses, inconsistent with the emission line 
fits. While for Kerr black holeaccretion disk models the black hole mass required to fit a given 
continuum depends sensitively on the inclination, the dependence in the Schwarzschild case is 
much less steep. For a Schwarzschild black hole one may therefore soften the constraint of small 
inclination, taking into account the intrinsic uncertainties in the mass determination and in the 
accretion disk -fitting method. This argument and the constraint of low black hole masses seem 
to infer that black holes in AGN (at least in the Seyfert 1 galaxies in our sample) tend to have 
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little or no rotation, being near-Schwarzschild . This subject will be elaborated in a separate work. 

8.3. Selection effects 

A major question that comes up when trying to apply the calibration result obtained from 
this sample to other AGNs without is whether the sample is representative. The objects in our 
sample - AGNs with adequate reverberation data - were selected by the criteria that make them 
"promising" for reverberation mapping: evidence or good chances to observe continuum and 
emission-linevariability. In most cases this meant some previous indication that the emission-lines 
varied. The evident selection criterion is therefore variability, in particular emission-line variability. 
However, variability is a feature common to many if not most AGNs, and it is difficult to identify 
an obvious way how such selection criteria can bias our results. If, for example, the more variable 
sources had systematically broader mean Hg line profiles, then the black hole masses derived 
using the coefficient from our sample would be underestimated. 

However, the line width enters in the same way in both mass estimation methods, only it is 
the mean profile in the photoionization method and the rms profile in the reverberation one, so 
it is the difference between the mean and rms profile width which should be considered, and that 
one shows no trend with luminosity (Fig. 6). We also note that the FWHM difference shows no 
preferred direction: it is positive for some objects and negative for others. 

Another question that may be posed is wheather the selection on objects on the basis of 
emission-line-variability introduces a bias at the parameter-space coverage: wheather our sample 
is representative of the AGN population. We have plotted our sample together with large samples 
of PG quasars and Seyfert 1 galaxies (both, narrow and broad line Seyferts) in the luminosity 
vs. H/3 line- width plane (Fig. 7). We can see that our sample provides a fair representation of 
the local AGN population, except of the high luminosity end. This limitation can be expected, 
as the BLR size is correlated with luminosity, so very luminous objects are likely to have larger 
BLRs, and hence longer response times, hence they would be more difficult for reverberation 
measurement. However, as Fig. 7 reveals no obvious correlation between luminosity and line 
width, this restiction is not likely to bias our results. 

8.4. Uncertainties of the reverberation method 

Although the BLR sizes determined from the lag are much more reliable than those given by 
the photoionization method, which involves much larger intrinsic uncertainties (see below), it has 
uncertain factors depending on the geometry of the emitting gas: 

• What is the error introduced by using moments of the transfer function rather than the 
full function? The shape of the cross-correlation function depends not only on the transfer 
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function, but also on the form of the continuum variations, so one can't read much from 

it. However, its first moment (the centroid) is rather less dependent on the form of the 
continuum variations. To REALLY address the problem correctly, one needs the transfer 
functions, but accurate transfer recovery is somewhat beyond the scope of what can be done 
with the data we have at present. 

• How is the lag time related to the geometrical distance of the line emitting gas from the 
illuminating continuum, and how does it depend on the geometry, in particular of a radially 
extended distribution? Also here there is a dependence on the continuum variations. For 

example, if the emitting region has an extended range of radii, Rin Rout the outer parts 

will not be able to respond to continuum changes with a timescale sorter than Rout/c, 
while the inner parts will. In that case a more adequate measure of the size may be obtained 
using the peak rater than the centroid of the cross-correlation function. 

More uncertainties are introduced by combining the lag with the line profile to give a mass 
estimate: 

• What measure of the line profile gives the appropriate velocity dispersion to combine with 
the lag? This becomes more complicated when the line emitting gas has a radially extended 
geometry, as the centroid lag, which is the responsivity weighted delay, may be differently 
weighted than the FWHM, which is the emissivity weighted line-of-sight velocity dispersion. 

• Does the width of line profile measure the velocity dispersion in the same gas the distance of 
which is measured by the lag? Using the rms rather than the mean profile reduces, but does 
not eliminate this uncertainty. 

• What is the relationship between the emissivity-wcightcd line-of-sight velocity measured by 
the line width and the central mass? This would depend on the shape of the orbits, which 
may introduce a systematic variation of a factor of a few. 

• What happens if the illuminating radiation does not originate directly from the central 
mass? In that case the lag and the line width would refer to different parts, the combination 
of which to give a virial mass estimate would depend on the specific geometry. 

All the above uncertainties are systematic, and introduce an error factor of a few, significantly 
larger than the random errors quoted in table 1. However, being systematic uncertainties, they will 
apply to all objects alike, altering the calibration coefficient but not effecting the linear correlation 
between the masses derived by the two methods. However, if there were luminosity-dependent 
effects, for example if the BLR geometry would change with L, or if the relation between the rms 
and the mean line width depended on L - of on M - this would change the slope of the correlation. 
As seen from Fig 6, such a systematic difference between the rms and mean FWHM in our sample 
is not present. 
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8.5. Uncertainties of the photoionization method 

The basic photoionization method apphed above to check the feasibihty of the cahbration 
makes several simphfying assumptions, which need to be examined in more detail or replaced by a 
more elaborate treatment. 

• Potential ambiguities arise as a result of the radial extent and ionization stratification of the 
BLR. The effect of these complications may be studied within the framework of the locally 
maximally emitting cloud model (Baldwin et al. 1995; Korista et al. 1995). 

• The physical parameters used to calculate the distance (density and ionization parameter ) 
may vary from object to object and differ also for different lines in the same object. The 
adjustable factor in the basic photoionization method depends on the values used for the 
"average" ionization parameter, density, and L/Lion- It is possible to use spectral data on 
different emission lines, to estimate the values of the ionization parameter and density for 
individual objects. 

• The mapping from the line width to the actual virial velocity depends on the distribution of 
the line emission in physical and velocity space, and how this is related to the line-width 
measure (e.g., FWHM, FWZI, intermediate width, skewness). The convolution of extended- 
emission BLR models with various possible dynamics of the BLR gas may be used to find 
an improved relation between the line profile and the effective velocity dispersion in the line 
emitting gas, induced by the central mass. 

• In extended BLR geometries, also the physical characteristics (such as density, differential 
covering area, emissivity) are likely to vary with radius, as indicated e.g. by two-zone 
photoionization modeling of several line ratios, which suggests a stratified BLR possibly 
extending into an "intermediate line region" (Brotherton et al. 1994). 

Compensating for these potential ambiguities in the photoionization method may give an improved 
photoionization method of virial mass determination, and may also reduce the systematic error 
sources in the photoionization mass estimates and the scatter in the Mphvs. Mrev diagram. 

Comparing masses and BLR sizes obtained by improved photoionization models to the 
reverberation - RMS values may in turn be used to test the photoionization BLR models. 

9. Summary 

We have estimated the virial mass of 19 AGN from the BLR size and emission-line width, 
using two different methods: the accurate reverberation mapping combined with the variable part 
of the emission-line, and the photoionization method, combined with the plain line width. We 
have shown that the two mass estimates for all objects in our sample are highly correlated, and 
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consistent with a linear relation, hence we were able to calibrate the photoionization method of 

mass and BLR size estimate. This should make possible mass estimates for virtually all AGN 
with H/3 line width measurements, which will be presented in a subsequent work. Our sample 
significantly enhances the available AGN central mass estimates: it has twice the number of 
objects previously available, it has better mass values, using the rms rather than the mean line 
profile, and all objects are uniformly treated. The calibration also enables us to estimate the 
ionizing luminosity, which cannot be directly observed. For our sample we find it to be on average 
10 times the monochromatic luminosity at 5100A. Combining the ionizing luminosities with the 
reverberation mass estimates our sample gives a mass-luminosity relation of M oc L"'®, with an 
Eddington ratio distributed in the range 0.01-0.3. 
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the hospitality of the Astronomy Department at UCLA. We are grateful for support of this work 
by the National Science Foundation through grant AST-9420080 to The Ohio State University. 
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Fig. 1. — The virial mass of NGC 5548 vs. FWHM-rms for H/? (6 observing seasons in the years 
1989-1996) and five other hnes in the years 1989 and 1993 (see text). Data taken from table 1 of 
Peterson and Wandel (1999). 

Fig. 2. — Virial mass obtained from the reverberation-rms data vs. the virial mass estimated by 
the photoionization method. The dotted line shows the best hnear {Mj-ev on Mph) fit weighted by 
the uncertainty in M^ev for each object. 

Fig. 3. — a. Reverberation centroid lags vs monochromatic luminosity at 5100A. b. Reverberation 
radius vs. BLR size estimated with the photoionization method. The dotted hne shows the best 
linear (r oc L-^/^ or Rrev oc Rph) fit weighted by the uncertainties in the lag. Formally, the 
uncertainties in the photoionization estimates are the uncertainties in L(5100A) which are very 
small. The intrinsic uncertainties in the photoionization estimates are not shown. 

Fig. 4. — Reverberation-RMS mass vs. monochromatic luminosity. The dotted line shows the best 
linear {Mrev o<i L) fit weighted by the uncertainties in M^ev 

Fig. 5. — Ionizing luminosity (derived from the lag, see text) vs. reverberation-RMS mass. The 
diagonal lines correspond to an Eddington ratio of 1 and 0.01. 

Fig. 6. — The fractional difference between the rms and mean FWHM(H/?) line profiles (FWHM 
difference divided by FWHM (rms) vs. luminosity. 

Fig. 7. — The reverberation AGN sample (circles) in the visual luminosity vs. FWHM(H/?) plane 
compared to other AGN: PG quasars (diamonds; from Boroson and Green, 1994) and Seyfert 1 
galaxies (triangles; from Boiler, Brandt and Fink, 1996). Open circles indicate the six years of data 
for NGC 5548. 



Table 1: Continuum Luminosity, FWHM H/? (mean and rms), and Continuum-Line Lags 



Name 


log ALa(5100A) 


WFWHM(mean) 


■yFWHM(rms) 


lag Tcent 




(erg s ^) 


/I — 1\ 
(km s ) 


(km s ) 


(days) 


3C 120 


43.59 


1910 


2210 ± 120 


43.8^^^:^ 


3C 390.3 


43.53 


10000 


10500 ± 800 


24.2l|^ 


Akn 120 


43.89 


5800 


5850 ± 480 


38.6j:^:f 


F9 


44.01 


5780 


5900 ± 650 




IC 4329A 


AC\ or* 

42.86 


5050 


5960 ± 2070 


1-4-2.9 


Mrk 79 


43.36 


4470 


6280 ± 850 


18 1+4-9 


MrkllO 


43.31 


1430 


1670 ± 120 


19.5l|^ 


Mrk 335 


43.53 


1620 


1260 ± 120 


i6.8ili 


Mrk 509 


43.91 


2270 


2860 ± 120 


79.3l|i 


Mrk 590 


43.44 


2470 


2170 ± 120 


2o.5t^:^ 


Mrk 817 


43.46 


4490 


4010 ± 180 


15.5lfi 


NGC3227 


42.07 


4920 


5530 ± 490 


10.9l5o'9 


NGC 3783 


42.79 


3790 


4100 ± 1160 


4.511? 


NGC4051 


41.47 


1170 


1230 ± 60 


6.5l«? 


NGC 4151 


42.62 


5910 


5230 ± 920 


3.0l{i 


NGC 5548^^ 


43.21 


6300 


5500 ± 400 


21.61^7 


NGC 7469^ 


43.54 


3000 


3220 ± 1580 


5-ol?:? 


PC 0804+762 


44.52 


3090 


3870 ± 110 


100.0lj|3 


PC 0953+414 


44.82 


2890 


3140 ± 350 


107.11^^:^ 



"For NGC 5548, the lag and the FWHM(rms) are derived from the whole data set of the years 1989 1996. 
''For NGC 7469, the cross-correlation lags are relative to the ultraviolet continnuum variations at 1315 A on account of 
the wavelength-dependent continuum time delays (see Wanders et al. 1997; Collier et al. 1998; Peterson et al. 1998b). 



Table 2: Reverberation BLR Sizes and Central Masses, Compared with Photoionization Sizes and 
Masses. The last two columns give the ionizing luminosity (derived from the lag, see text) and the 
corresponding Eddington ratio 



Name log log lag log Mph log M^et; Mrev{W^ Mq) log Lion log{Lion/LEdd) 



3C120 


0.92 


1.64 


6.86 


7.49 


q 1+2.0 


45.03 


-0.57 


3C 390.3 


0.89 


1.38 


8.26 


8.59 


39.11}^ 


44.51 


-2.19 


Akn 120 


1.07 


1.59 


7.97 


8.29 




44.92 


-1.48 


F9 


1.13 


1.23 


8.03 


7.94 


8.7tli 


44.21 


-1.84 


IC4329A 


< 0.56 


0.15 


7.34 


< 6.86 


< 0.73 


< 42.04 


< -2.93 


Mrk79 


0.81 


1.26 


7.48 


8.02 


10.5+f° 


44.26 


-1.87 


MrkllO 


0.78 


1.29 


6.46 


6.91 


0.80l°i^ 


44.33 


-0.69 


Mrk335 


0.89 


1.23 


6.68 


6.58 


0.39t°it 


44.20 


-0.49 


Mrk509 


1.08 


1.90 


7.17 


7.98 




45.54 


-0.54 


Mrk590 


0.85 


1.31 


7.00 


7.15 




44.37 


-0.89 


Mrk817 


0.86 


1.19 


7.53 


7.56 




44.13 


-1.54 


NGC 3227 


0.16 


1.04 


6.92 


7.69 




43.82 


-1.98 


NGC 3783 


0.52 


0.65 


7.05 


7.04 


1 1+1-1 

-■^•-L-i.o 


43.05 


-2.10 


NGC 4051 


-0.14 


0.81 


5.37 


6.15 


"-'•-'-^-0.09 


41.57 


-0.84 


NGC 4151 


0.44 


0.48 


7.35 


7.08 


1 9+0.8 
-■-•^-0.7 


42.70 


-2.49 


NGC 5548 


0.73 


1.26 


7.70 


7.83 


6.8li| 


44.27 


-1.83 


NGC 7469 


0.90 


0.70 


6.87 


6.88 




43.14 


-1.86 


PG 0804+762 


1.39 


2.00 


7.74 


8.34 


21.9+^1 


45.75 


-0.70 


PG 0953+414 


1.54 


2.03 


7.83 


8.19 




45.81 


-0.49 



Table 3: Reduced and probabilities (in parenthesis) for the three fitting schemes, constant mass 
(columns 3-4), lagoc {FWHM)'^ (y=-2x-|-a, columns 5-6), and the best fit (y=-bx-|-a, columns 7-9). 
The first column indicates the data set, the second one (n) - the number of points. 



Data set 


n 


M7 


xHm) 


a 




b 


a 


X^{—bx + a 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


All lines, RP 89-96 


14 


6.1 


1.09(0.35) 


8.85 


2.10(0.04) 


-2.05±0.13 


8.85± 0.49 


1.78(0.05) 


All lines, H/3 89 


9 


5.5 


0.92(0.49) 


8.58 


1.57(0.15) 


-1.98±0.13 


8.50± 0.51 


1.22(0.29) 


All lines, H/3 average 


9 


6.0 


0.87(0.51) 


8.62 


1.44(0.19) 


-2.22±0.18 


9.51± 0.69 


0.96(0.45) 


H/3 89-96, centroid 


6 


6.6 


1.72(0.15) 


8.66 


3.51(0.01) 


-0.76±0.39 


4.06± 1.43 


0.93(0.45) 


H/? 89-96, peak 


6 


6.6 


0.58(0.68) 


8.66 


1.67(0.15) 


-1.47±0.27 


6.67± 1.00 


0.66(0.62) 



NGC 5548 - various lines 
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